AMENDMENTS TO THE CLAIMS 

This listing of claims will replace all prior versions, and listings, of claims 
in the application: 

Listing of Claims; 

1 1 . (Currently amended) A method for bounding the solution set of a 

2 system of linear equations Ax = b, wherein A is an interval matrix and b is an 

3 interval vector, the method comprising: 

4 receiving the system of linear equations Ax = b; 

5 storing Ax = b in a memory in a computer system; 

6 preconditioning the set of linear equations Ax = b by multiplying tlii-ough 

7 both side of the linear equations b y a matrix B to produce a preconditioned set of 

8 linear equations BAx=Bb, wherein the set of linear equations is a representation 

9 of a global optimization problem; 

1 0 substituting Mo = B A and r = Bb to produce Mqx = r; 

1 1 widening the matrix Mo to produce a widened matrix M, wherein the 

12 midpoints of the interval elements of M form the identity matrix; and 

13 using M and r to compute the hull a hull h of the system Mx = r, which 

14 boimds the solution set of the system Mox = r; 

15 wherein the interval operations involved in bounding the solution set are 

16 perfonned using a special-purpose interval arithmetic unit configured to perform 

17 interval arithmetic operations . 

1 2. (Previously presented) The method of claim 1 , wherein the method 

2 further comprises computing the matrix B by: 
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3 computing an approximate center Ac of the interval elements of matrix A; 

4 and 

5 forming B by computing an approximate inverse of Ac, B = (Ac)'^ 
1 3 (Canceled). 

1 4. (Previously presented) The method of claim 1, further comprising 

2 assuring that sup(ri) ^ 0 by changing the sign of r/ [[(]]and X/[[)]] if necessary, 

3 v^herein r, is an element of r. 

1 5. (Original) The method of claim 1 , further comprising: 

2 determining if M is regular; and 

3 using the Gauss-Seidel process for computing the hull h if M is not 

4 regular. 

1 6. (Currently amended d) A computer-readable storage medium storing 

2 instructions that when executed by a computer cause the computer to perform a 

3 method for bounding the solution set of a system of linear equations Ax = b, 

4 wherein A is an interval matrix and b is an interval vector, the method 

5 comprising: 

6 receiving the system of linear equations Ax = b; 

7 storing Ax = b in a memory in a computer system; 

8 preconditioning the set of linear equations Ax = b by multiplying through 

9 both side of the linear equations by a matrix B to produce a preconditioned set of 

1 0 linear equations BAx=Bb, wherein the set of linear equations is a representation 

11 of a global optimization problem; 

12 substituting Mo = BA and r = Bb to produce Mqx = r; 
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widening the matrix Mo to produce a widened matrix M, wherein the 
midpoints of the interval elements of M form the identity matrix; and 

using M and r to compute the hull a hull h of the system Mx = r, which 
bounds the solution set of the system Mox = r; 

wherein the interval operations involved in bounding the solution set are 
performed using a special-purpose interval arithmetic unit configured to perform 
interval arithmetic operations . 

7. (Previously presented) The computer-readable storage medium of claim 
6, wherein the method further comprises computing the matrix B by: 

computing an approximate center Ac of the interval elements of matrix A; 

and 

forming B by computing an approximate inverse of Ac, B = (Ac)"^ 

8. (Previously presented) The computer-readable storage medixim of claim 
6, wherein using M and r to compute the hull h involves: 

forming P as an inverse of the left endpoint of M; 
forming c, = 1/(2P// - 1) for / = 1, . n; 

forming z, = {in/{ri)-\-sup(n))Pii - elVsup{r\ wherein el'i^ a unit vector in 
which the /-th element is 1 and other elements are 0, and wherein r, is an element 
of r; 

setting inf(hi) = c^i if z, > 0; 
setting inf{hi) = Zi if z, ^ 0; and 
setting sup(h) = Fsup{r), 

9. (Previously presented) The computer-readable storage medium of claim 
6, wherein the method further comprises assuring that sup(n) ^ 0 by changing the 
sign of r/ [[(]]and x/[[)]] if necessary, wherein r, is an element of r. 
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1 10. (Original) The computer-readable storage medixim of claim 6, wherein 

2 the method fiirther comprises: 

3 determining if M is regular; and 

4 using the Gauss-Seidel process for computing the hull h if M is not 

5 regular. 

1 11. (Currently amended) An apparatus that bounds the solution set of a 

2 system of linear equations Ax = b, wherein A is an interval matrix and b is an 

3 interval vector, comprising: 

4 a receiving mechanism configured to receive the system of linear 

5 equations Ax = b; 

6 a special purpose arithmetic unit configured to perform intei-val i 

7 computations; 

8 a storing mechanism configured to store Ax = b in a memory in a ! 

9 computer system; ^ 

1 0 a preconditioning mechanism within the special purpose arithmetic unit ^ 

1 1 that is configured to precondition the set of linear equations Ax = b by 

1 2 muhiplying through both side of the linear equations by a matrix B to produce a { 

1 3 preconditioned set of linear equations B Ax=Bb, wherein the set of linear 

14 equations is a representation of a global optimization problem; 

1 5 a^substituting mechanism within the special purpose arithmetic unit tliat is 

1 6 configured to substitute Mo = BA and r = Bb to produce Mqx = r; 

1 7 a widening mechanism within the special purpose arithmetic unit that is 

1 8 configured to widen the matrix Mo to produce a widened matrix M , wherein the 

1 9 midpoints of the interval elements of M form the identity matrix; and 

20 a hull computing mechanism within the special pui-pose arithmetic unit 

21 that is configured to use M and r to compute tbe-faril -a hull h of the system 

22 Mx = r, which bounds the solution set of the system Mqx = r; 
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wherein the interval operations involved in bounding the solution set are 
performed using the special-purpose interval arithmetic unit . 



1 12. (Previously presented) The apparatus of claim 1 1 , wherein the 

2 preconditioning mechanism is configured to: 

3 compute an approximate center Ac of the interval elements of niatrix A; 

4 and to 

5 form B by computing an approximate inverse of Ac, B = (Ac)'^ . 

1 13. (Previously presented) The apparatus of claim 1 1 , wherein the hull 

2 computing mechanism is configured to: 

3 form P as an inverse of the left endpoint of M; 

4 form d = \l{2Pii - 1) for / = 1, . . n\ 

5 form Zi = {infiri)+sup{ri))Pii - elVsup{Y), wherein elis a unit vector in 

6 which the /-th element is 1 and other elements are 0, and wherein r, is an element 

7 ofr; 

8 set infih^) = C/Z, if z, > 0; 

9 set infih^ = z, if z/ < 0; and to 
1 0 set sup(h) = Vsup{Y). 

1 14. (Previously presented) The apparatus of claim 1 1, wherein the 

2 preconditioning mechanism is configured to assure that sup{r^ > 0 by changing 

3 the sign ofr, [[(]]and X/[[)]] if necessary, wherein is an element of r. 

1 15. (Original) The apparatus of claim 1 1, wherein the preconditioning 

2 mechanism is configured to: 

3 determine if M is regular; and to 

4 terminate the process of computing the hull h if M is not regular. 
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16. (Currently amended) A method for bounding the solution set of a 
system of linear equations Ax = b by multiplying through both side of the linear 
equations by the matrix B to produce a preconditioned set of linear equations 
BAx=Bb, wherein the set of linear equations is a representation of a global 
optimization problem, the method comprising: 

receiving the system of linear equations Ax = b; 
storing Ax = b in a memory in a computer system; 
substituting Mo = BA and r = Bb producing Mqx = r; 
assuring that sup(r^ ^ 0 by changing the sign of r, (and x,) if necessary; 
widening the matrix Mo to produce a widened matrix M, wherein the 
midpoints of the interval elements of M form the identity matrix; and 

using M and r to compute the hull a hull h of the system Mx = r, which 
bounds the solution set of the system Mox = r by, 

forming P as an inverse of the left endpoint of M, 
forming c, = 1/(2P„ - 1) for / = 1, . . 
forming z, = {infiri)-^sup(ri))Pii - el^supir), 
wherein el\^ a unit vector in which the /-th element is 1 
and other elements are 0, and wherein is an element of r, 
setting infihi) = CiZj if z, > 0, 
setting iri/{hj) = z, if z, < 0, and 
setting sup(h) = Vsup(r)i 
wherein the interval operations involved in bounding the solution set are 
performed using a special-purpose interval arithmetic unit configured to perform 
interx^al arithmetic operations . 

17. (Original) The method of claim 16, further comprising: 
determining if M is regular; and 
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3 using the Gauss-Seidel process for computing the hull h if M is not 

4 regular. 

1 18. (Previously presented) The method of claim 1 6, wherein the method 

2 further comprises computing the matrix B by: 

3 computing an approximate center Ac of the interval elements of matrix A; 

4 and 

5 forming B by computing an approximate inverse of Ac, B = (Ac)'* . 
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